Canted ground state in artificial molecules at high magnetic fields. 
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We analyze the transitions that a magnetic field provokes in the ground state of an artificial 
homonuclear diatomic molecule. For that purpose, we have performed numerical diagonalizations 
for a double quantum dot around the regime of filling factor 2. We present phase diagrams in terms of 
tunneling and Zeeman couplings, and confinement strength. We identify a series of transitions from 
ferromagnetic to symmetric states through a set of canted states with antiferromagnetic couping 
between the two quantum dots. 



PACS numbers: 73.40.Hm 

It is well known Q, that the ground state (GS) of 
homonuclear diatomic molecules, as C2, O2 or N 2 , suf- 
fers a series of transitions when the distance between the 
nuclea is varied. Particularly interesting are the singlet- 
triplet transitions, which determine the magnetic prop- 
erties of gases formed with those molecules. Since the in- 
tramolecule interatomic distance is, essentially, constant 
it is not easy to control the physical mechanisms respon- 
sible for such a behavior. Recently, double quantum dots 
(DQD), which can be thought of as artificial homonu- 
clear molecules (|||, have opened up new possibilities, 
that should allow both a better understanding of the 
physics of these systems and the tailoring of their mag- 
netic properties, as some parameters, which are fixed in 
natural molecules, can be continuously varied in DQD. 
For instance, it is possible to change potential barrier 
heights, thus varying the tunneling rate between the two 
artificial atoms, without altering the interdot distance, 
which determines the electron-electron interaction. Also, 
the interdot distance is a parameter that can be exter- 
nally controlled by appropriate design of the nanostruc- 
ture defining the DQD. 

The use of an external magnetic field B to produce 
transitions in the GS is a particularly interesting tool to 
study artificial molecules: apart from the effect of Zee- 
man energy, the magnetic field introduces an additional 
length scale (the magnetic length) that, even for mod- 
erate fields, can be of the order of the molecular dimen- 
sions in these artificial systems. In this work we analyze 
the symmetric-ferromagnetic transition which occurs in 
DQD's with N electrons as the magnetic field is increased. 
Some previous works |^,^| have dealt with DQD in mag- 
netic fields, but in a totally different regime: there, the 
spin degree of freedom was not included, i.e., those works 
only considered magnetic fields high enough for the sys- 
tem to be always in the ferromagnetic phase. 



We restrict ourselves to homopolar molecules, i.e., 
DQD with N even where the electrons are three dimen- 
sionally confined: in the z direction they can reside in 
either one of two layers (top (T) or bottom(B)), sep- 
arated by a distance d, while in the x — y plane elec- 
trons are confined by the presence of a two-dimensional 
parabolic potential V(r) = m*uj r 2 /2 (m* is the elec- 
tron mass, r is the in-plane distance) and to a range 
of magnetic fields such that the DQD has a filling fac- 
tor close to v = 2. This last restriction is motivated 
by the symmetric-ferromagnetic transitions found in in- 
finite multicomponcnt Quantum Hall (QH) systems close 
to filling factor v = 2 @. 

We work in the symmetric gauge, with a set of states 
restricted to the lowest orbital Landau level, i.e. single 
particle wavefunctions without nodes in the radial direc- 
tion. In this basis set, the Hamiltonian of the DQD is 

H=aM- A Z S Z - ( c L*t c ™°b + c^ CTS c mtrT ) 
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where co c = eB/m* and a = h([uj 2 + -iuj 2 ] 1 / 2 - lu c )/2. 
m and M are the third components of single particle 
and total angular momentum respectively, a and S z 
are the third components of single particle and total 
spins respectively and A a layer index, equal to T or 
B. A z — gfJ-sB is the Zeeman coupling with g being 
the Lande g-factor and \xb the Bohr magneton. A sas is 
the single particle energy gap between symmetric (s) and 
antisymmetric (as) combinations of T and B single par- 
ticle states. V AA =e 2 /(er) and V AA ' =e 2 /(eVr 2 + d?) for 
A ^ A' are Coulomb interaction potentials, e being the 
electron charge and e the dielectric constant. In all our 
results, the energies are given in units of e 2 /(els) with 
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Ib = [Ti/{m*[oj 2 c +<iuj 2 ] 1 / 2 )] 1 / 2 being the magnetic length. 
The eigenstates of the DQD are a function d, A z , A sas 
and a. The spectrum is separated in subspaces labeled 
by the quantum numbers (M, S z , Pi) where Pj = (— 1) ///2 
is the parity of the isospin I = N s — N as with N s , N as 
being the number of electrons in symmetric and anti- 
symmetric states respectively. It should be noted that, 
for d 7^ 0, intcrdot and intradot Coulomb interactions are 
different. In this case, electron-electron interaction mixes 
states with different isospin, and I is not a good quan- 
tum number. However, the parity Pi is always a good 
quantum number because it is related to a symmetry op- 
eration of the problem: the reflection 7Z with respect to 
the midplane between the two layers. The application 
of 1Z to any state (not necessarily an eigenstate) with a 
well defined isospin would result in the same state mul- 
tiplied by (-1)^" = (-l)( Ar ~ / )/2 which, apart from the 
unessential constant (— l) N ^ 2 , is precisely Pj times the 
state. This shows that 1Z applied to any state without 
a well defined isospin would still produce the same state 
(except for a sign), provided the initial state had weight 
in different isospin subspaces with the same isospin par- 
ity. There is another symmetry operation in the problem: 
the inversion with respect to the midpoint between the 
centers of the two quantum dots (QD). We do not pay 
special attention to it in the discussion of our results, 
because it does not add any information of relevance. 

We have diagonalized numerically the Hamiltonian (|l|) 
for N — 6 and N = 8. As similar structure is found in 
both cases, only results for N = 8 are presented here. 

The analysis of the results will be helped by the large 
experience accumulated in the knowledge of the GS of 
a multicomponent QH systems Q. When electrons are 
confined in a double layer (DL) in the regime of global 
filling factor v = 2, they have as degrees of freedom a 
(=T, I) and layer index A. Those degrees of freedom pro- 
voke a rich phase diagram in terms of tunneling and Zee- 
man couplings. Obviously, when Zecman splitting A z is 
much larger than tunneling splitting A sas , the system 
prefers a ferromagnetic GS \F) in which electrons oc- 
cupy all the symmetric and antisymmetric states with 
spin \. In the opposite limit A SQS >> A z , the GS 
\S) corresponds to electrons fully occupying the sym- 
metric states with both f and J. spins. The character 
of the transition between these two extreme situations is 
controlled by electron-electron interaction effects in the 
four-dimensional space of degrees of freedom. A very 
attractive proposal was made for a GS, labeled as 
canted, which connects continuously between the \F) and 
\S) limiting cases. The canted state is a ferromagnet in 
the field direction while, for the direction perpendicular 
to the magnetic field, spins in different layers have anti- 
ferromagnetic correlations. The properties of this state 
has been studied by a microscopic Hartree-Fock theory, 
a long wavelength field theory based on the quantum 
0(3) nonlinear sigma model and a bosonic spin theory 



Moreover, exact numerical diagonalizations for a 
small number of electrons in a DL with spherical shape 
]l5|] show that, in a translationally invariant system, the 
canted phase survive quantum fluctuations; although the 
domain, in the (A z , A$as) parameter space, in which it 
is the GS is much narrower than what Hartree-Fock ap- 
proximation predicts. From the experimental side, such 
a phase is consistent with the available information in 
DL by inelastic light scattering p6[, magnetoresistance 
]l7t and capacitance spectroscopy |18[|. 

In our calculations in DQD, we also find that, for large 
A z and moderate A sas , the DQD has a IF) GS while for 
large A sas and moderate A z the GS is \S). As in the infi- 
nite system, the DQD presents a transition between these 
two regimes through a set of more complicated states. 

In order to identify a GS obtained from numerical diag- 
onalizations as a canted state, a difficulty arises from the 
fact that the eigenstates of the DQD have a well defined 
third component S z of the spin while the mean field wave 
function \C MF ) for the canted state is not an eigenstate 
of S z . However, we can restore the broken symmetry 
of the mean field states by projecting on subspaces with 
well defined S z = N/2 - n: 
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where n is an integer number. In Eq. (|2j), (f> is the 
angle defining a particular canted state \C MF ((/))) = 
(<f>i_ (0), $2- (0)) in which, for filling factor 2, the elec- 
trons are occupying two type of states 
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The notation is the one used by Das Sarma et al. ^ with 
the states (||) given in the basis of (s, t), (as, I), (s,|), 
(as, "f). The angles 6i depend on the Hamiltonian param- 
eters being B\ = 2 = for the ferromagnetic state \F) 
and 0i = and 2 = 7r for the symmetric state \S). After 
some algebra, one gets the projections (||) of the canted 
state written as 



\cif F ) = C n 
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\F) (4) 



C n being a normalization constant. The coefficients in 
Eq.(^|) are m-independent in the infinite system because 
all single particle states are degenerate in the lowest Lan- 
dau level. However, in a parabolic QD, translational in- 
variance is broken, the single particle energy depends on 
m and the corresponding canted eigenstates of the DQD 
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should have a structure similar to Eq.(|4|) although the 
coefficients 9™ and 6™ could have a dependence on m. 

Fig. [l] contains our results in a phase diagram showing 
the quantum numbers (M, S z , Pi) of the GS for d = Ib 
and N = 8. In order to work in the regime v — 2, we 
chose ol = 0.2e 2 /{els)- In the region to the left corre- 
sponding to small A sas , the GS, labeled as (12,4,1) is a 
ferromagnetic state given by a single Slater determinant 

\F) = nm=o c L s ,T c La S ,r There is another large fer- 
romagnetic region labeled as (13,4,-1) with a GS which 
is practically (more than 98%) one Slater determinant 

i^') = nt=o c L,^Tnm=o c La S ,r This mcans that for 

an increasing tunneling, the ferromagnetic state \F) suf- 
fers an edge reconstruction by a charge instability which 
preserves the ferromagnetic character giving \F'). In 
the lower region corresponding to small A z , the GS, la- 
beled as (12,0,1), is a symmetric state given by prac- 
tically (more than 98%) one Slater determinant \S) — 
rim=o c m s T c m s i- When the Zeeman coupling increases, 
there is a new symmetric GS, labeled as (13,1,1), which is 

also practically given by \S') = nl =0 cL, s , T Xlm=a c L,s,l- 
The increase of A z provokes an edge reconstruction which 
involves a charge-spin excitation giving a new GS which 
essentially preserves the symmetric character. 

The most interesting part of the phase diagram cor- 
responds to the narrow regions separating ferro-like re- 
gions from symmetric-like regions. The regions (12,3,-1), 
(12,2,1) and (12,1,-1) contain canted states \C n ) (corre- 
sponding to n = 1,2,3 respectively) while the regions 
(13,3,1) and (13,2,-1) contain canted states \C' n ) (corre- 
sponding to n = 1,2 respectively). This identification is 
made through the wavefunctions which turn out to have 
the functional form given by Eq. (^) with coefficients 9™ 
and 0™ almost independent on m at the center of each 
QD. Since these regions are well described by states (^), 
they must be understood as the projection to well defined 
quantum numbers of canted states with an antiferromag- 
netic tilting of the spins. 

Apart from the states shown in Fig. [j], some other GS, 
in particular with different edge reconstructions, could 
appear when the parameters are varied in the calcula- 
tion. In fact this is the case: as shown in figures |^ and 
|3|, showing phase diagrams for d — Ib, N — 8, when 
A z and a are varied for fixed A sas = O.le /(els) and 
A sas = 0.2e 2 / (els) respectively. In the central part of 
Fig. | the sequence (12,4,1), (12,3,-1), (12,2,1), (12,1,-1) 
and (12,0,1) corresponds to the already discussed tran- 
sition |F) — > \C n ) — > \S) with n = 1,2,3. Contiguously 
to its left, for smaller values of a, the sequence (13,4,-1), 
(13,3,1), (13,2,-1) and (13,1,1) corresponds to the tran- 
sition \F') -> \C' n ) -> \S f ) with n = 1,2. To the right 
of those regions, the increase in a implies a higher con- 
finement and, consequently, a larger concentration of the 
electrons in the center of the DQD which, in turns, im- 
plies a decrease of M. Edge reconstruction occurs in the 



leftmost region, i.e., smallest a, where a reduction of the 
confinement provokes the spreading of the electron den- 
sity away from the center. In Fig. || the tunneling is so 
large that the sequence (12,4,1), (12,3,-1), (12,2,1) and 
(12,1,-1) does not appear and only the symmetric state 
1 5) = (12,0,1) is observed. Apart from this, the whole 
structure is similar to that shown in Fig. |[ 

In summary, we have analyzed the transitions that a 
magnetic field provokes in the GS of an artificial homonu- 
clear diatomic molecule. For that purpose, we have car- 
ried out numerical diagonalizations for a DQD for filling 
factors close to v = 2. The resulting phase diagrams, 
are understood to the light of previous experience on 
the GS of an infinite DL at the same regime. When 
the different parameters are varied, a series of transi- 
tions \F) — ► \C n ) — > | S) from ferromagnetic to symmet- 
ric states are identified, through a set of projections (into 
subspaces with well defined quantum numbers) of canted 
states. Such canted states have ferromagnetic correla- 
tions in parallel (to magnetic field) direction and antifer- 
romagnetic correlations in the perpendicular components 
of spins in different QD. The electron-electron correlation 
that leads to the canted states in infinite double layers for 
v = 2 survives, in this regime, to the presence of edges 
and is present even when there is edge reconstruction. 
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FIG. 2. Phase diagram showing the subspace (M, S z , Pi) 
containing the GS for d = Ib, A sas = 0.1e 2 /(eZs) and N — 8. 
Energies a and A z are given in units of e 2 /(els)- 
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FIG. 1. Phase diagram showing the subspace (M,S z ,Pi) 
containing the GS for d — Ib, a = 0.2e 2 /{eIb) and N — 8. 
Energies A 3as and A z are given in units of e 2 /(eIb)- 
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FIG. 3. Phase diagram showing the subspace (M,S z ,Pi) 
containing the GS for d = Ib, A sas = 0.2e 2 /(eZs) and N = 8. 
Energies a and A z are given in units of e 2 /(els)- 
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